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In a recent study [Phys. Rev. E 60, 6530 (1999)], Trizac and Raimbault showed that the effective 
pair interaction between like charged colloids immersed in a cylindrically confined electrolyte remains 
repulsive even when the size of the micro-ions or the finite longitudinal extension of the confining 
cylinder are taken into account. Contrary to their claim, we argue that the case of finite longitudinal 
confinement doesn't always generate repulsive interactions and to illustrate this point we also provide 
a simple example. 

PACS numbers; 05.70.Np, 64.10+h, 82.60. Lf, 82.70.Dd 

There have been a number of recent theoretical efforts towards understanding the mechanism behind the long range 
attraction that is experimentally observed [l]-|^ when like charged colloids are immersed in a confined electrolyte. The 
fact that this interaction returns to a purely repulsive one as soon as the confining surface is removed , suggested 
that the attractive part could be obtained within the framework of the Poisson-Boltzmann theory by properly taking 
into account the influence of the confining surface. However, Neu [|| and Sader and Chan using a general model 
of cylindrical confinement, have rigorously proven that this is not the case. 

In their paper, Trizac and Raimbault |Q extended the aforementioned studies [||J^ to include the influence brought 
upon the sign of the interaction by the size of the micro-ions and the finite longitudinal extension of the confining 
cylinder (case hereafter referred to as complete confinement). Although we agree with their conclusion that the 
interaction remains repulsive when the specific size of the micro-ions is taken into account through a modified Poisson- 
Boltzmann equation, we don't agree with a similar conclusion for the case of complete confinement. In this comment, 
we show that this controversy is brought about by an inaccuracy in their calculation and that in general, colloids 
immersed in a completely confined electrolyte don't always repel. We also provide a simple example to illustrate this 
point. 

To be more specific, the issue has to do with the use of Green's theorem in two dimensions (2Z)). If is a closed 
domain in (for example in the xy plane) and 4>(x,y) and '4'(x,y) are scalar fields defined on R^, then Green's 
theorem in 2D is: 
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where n (which is contained in the xy plane) represents the outer unit normal to the boundary dW of W and V^y 
is the 2D gradient operator. 

With this in mind and making use of the same notation and assumptions as in , we obtain a different expression 
for Eqs. (14) and (15) in that paper, where Green's theorem was inaccurately applied. For clarity, let us first separate 
the z component of the gradient in the equation that generates the results (14) and (15) in 0: 
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where V is the 3-D gradient operator and tp possesses mirror symmetry with respect to the plane z = 0, as in [Q. 
The reader should remember that the notation Jp^y actually stands for the two integrals over the domains of the 
xy plane that are interior and exterior, respectively, to the confining cylinder. The dimensionality of these integrals 
requires the use of the 2D version of Green's theorem, unlike it was inaccurately done in R], where the 3-D version 



was employed instead. This issue is clarified if we consider the more explicit repr esentation shown in Eq. (0.4), where 



we can correctly apply the 2-D version of Green's theorem [as stated in Eq. (OA)] to its first term to get 
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where in going from Eq. (0.5) to Eq. (0.6) we used the dielectric boundary condition D'^ — = a (and the fact 
that Jg^^di {ipz=o — iPz=l) — i^dS Ez) and also formed the 3D Laplacian and employed Eq. (10) in 0. 

In fact, Eqs. (3.5;) and (0.6) are the correct forms for Eqs. (14) and (15), respectively, in m. After adding all the 
terms, the correct expression for the z component of the force (Eq. (16) in [Q) should read: 
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Although one can use the same argument as in to prove that the first two integrals in Eq. ( p.T] ) are always 
positive, the third one (which is different from that obtained by Trizac and Raimbault) can in fact be negative and 
sometimes overcome the first two. Before we show how this can happen, let us mention another equivalent expression 
for Fz, which can be obtained by recasting in a different manner the second integrand of Eq. (13) in [Q, namely 

(D-E)^^g — ( D-E )^_r = — e (E^^o — Ez=l)'^ — 2Eiz=q - {T)z=l — TOz=o) and then following the same steps used in 
deriving Eq. (0.7): 
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The reader can easily check that Eqs. ( p.7| ) and (3.6) are equivalent, by subtracting them from one another and 
then using Eq. (10) in |^], Green's theorem in 2D and the dielectric boundary condition to prove that the result is 
zero. 

We can now propose a simple example in order to prove that the complete confinement doesn't always generate 
repulsive interactions. As before, we assume that the particles are immersed in an electrolyte of permittivity e, confined 
by an infinite and uniformly charged cylindrical surface (of arbitrary cross section), while the medium outside the 
cylinder has a permittivity e' and also contains an electrolyte (most general case). The only requirement is that ip 
possesses mirror symmetry with respect to the plane z = 0. In our example the particles Si and S2 are taken to 
be flat surfaces with the same surface charge density a as the confining cylinder while the micro-ions are able to 
equilibrate at all times (by being in contact with a reservoir and possibly permeating through the colloidal particles). 
The walls A and B that realize the complete confinement are surrounded on both sides by electrolyte and also have 
the same uniform surface charge density a as the confining cylinder (this way the complete confinement is realized by 
a surface of constant a). To compute the z component of the force on the colloid S2, we follow the same procedure 
outlined in only that instead of choosing the surface to integrate over in Eq. (11) that defines Fz in Q], we 
choose the surface E^, midway between the colloid S2 and the wall B (see Fig. |l|). Keeping in mind the correction 
made in this comment, the result ends up to be identical to Eq. (^^), with the only exception that all the indices 
z = L are replaced hy z — d: 
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It is easy to see that the last two integrals in Eq. ( |0.9| ) approach zero as the separation 2L between the walls is 
increased to infinity and the distance 2(L — d) between each particle and the adjacent wall is kept constant. Indeed, in 
this situation, ip becomes symmetric with respect to the plane z = d (due to our choice for the shape of the colloidal 
particle) and therefore {Ez)^^^ — > 0, while at z = (infinitely far away from the particles and the walls), the problem 
of solving for -0 becomes essentially a 2D one (independent of z) and therefore d^ip/dz^ — > for any n > 1. At 
the same time, as argued in the first two integrals in Eq. ( |0.9| ) are always positive (the first integrand is positive 
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due to the global convexity of P{ip)) and moreover, reach some finite asymptotic values in this limit case. Since both 
integrals contribute with a minus sign in Eq (^^), it becomes clear that < when L — > oo and L — d is constant. 

Therefore, within the approximations of our example, we conclude that for any fixed distance 2{L — d) between the 
colloidal particle and the adjacent wall, there is a finite, maximum separation between the walls 2Lo, at which the 
particles are in equilibrium (F^ — 0) and beyond which F^ becomes negative and remains so as L increases. Although 
this is in contradiction with the result obtained earlier by Trizac and Raimbault |0] that F^ should always be positive 
under complete confinement, we have shown that the controversy stems from an inaccuracy in their calculation. It is 
therefore possible that completely confined colloidal particles don't always repel, mainly due to the stronger interaction 
with the walls (or conceivably other particles |T^]) that are responsible for the finite longitudinal confinement and 
break the translational symmetry of the confined space. Although this type of calculation is not suitable to yield an 
estimate on Lq, it would nevertheless be of high interest if this could be done through computer simulations. 
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FIG. 1. Longitudinal section along the confining cylinder showing the two particles Si and S2 and the walls A and B 
responsible for the complete confinement. The plane z = d is midway between the particle S2 and the wall B. 
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